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A small-signal (linearized) theory of discrete-charge-transfer-device 
performance is presented for the case of incomplete charge transfer. 
Specifically, the dispersion relation is derived which relates the charge- 
transfer efficiencies presently characterizing these discrete (in space 
and time) devices to the usual measures of device or transmission-line 
performance based on the attenuation , dispersion, phase velocity, etc., 
of sine waves. In a more general sense this emphasizes the applicability 
of conventional signal theory to these new devices. The impulse 
solution or Green's function is then shown to be the equivalent of a 
bivariate distribution in probability theory. More generally the utility 
of ( deter ministically interpreted) probability theory is emphasized 
by showing the equivalence of a general s?7iall- signal theory to a 
random-walk process. 

I. INTRODUCTION 

In an important new class of discrete-eharge-transfer devices includ- 
ing charge-coupled devices 1 (CCD's), bucket-brigade shift registers, 2-3 
and other shift-register or image-detection or display devices, ex- 
ternally applied time-dependent voltages step captive charge along 
a chain of equivalent discrete storage stations. In some of these 
devices the charge transfer is imperfect with a fraction of the charge 
failing to advance and a fraction lost altogether during each step. 
Explicit expressions are constructed here for the dispersion relations 
and Green's functions which describe this imperfect performance 
under conditions when the fractions of charge that go astray can be 
described by constant parameters characteristic of the particular 
device. 

The theory of analog signal processing is based on the properties 
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of sine waves, and analog devices are conventionally characterized 
by the dispersion, attenuation, etc., which they cause. A second 
standard method of characterization is based on the distortion and 
attenuation of pulses. In contrast, discrete-charge-transfer devices 
are presently characterized by their charge-transfer efficiencies 
(fractions). The solutions given here were chosen to facilitate the 
application of conventional signal processing theory to these new 
discrete devices. That is, they show the equivalence of these three 
methods of characterization in the .small-signal limit and provide the 
appropriate interrelating formulas. Specifically, Section III treats the 
sinusoidal representation, while Section IV considers pulses. 

Although the physical interpretation of our equations is determinis- 
tic, it is also an objective in Section IV to show the direct applicability 
of many established results of probability theory. In the Appendix our 
basic equation is re-derived as the simplest nontrivial example of a 
discrete small-signal theory. In turn, the small-signal theory is seen 
to be the deterministic limit of a random-walk process. 

Whether the mode of device operation be digital or analog, one is 
ordinarily interested in utilizing the full information capability; 
i.e., maximum bandwidth or wavelengths approaching twice the sta- 
tion separation in shortness. Thus, we at no time approximate the 
discrete equation by a (continuous) differential equation. 

In summary our objective is to emphasize by way of a case of 
present interest that standard methods of device characterization, as 
well as established probability theory, can be applied to discrete 
charge-transfer devices. 

II. BASIC EQUATION 

Discrete-charge-transfer devices are often, to a good approximation 
if not exactly, discrete in both space and time. That is, the informa- 
tion-bearing charge is moved in discrete bursts in time from one 
spatially discrete storage station (e.g., capacitor or potential well) to 
the next along a line of stations. Consider q x ,t the charge in station 
x at time t where x and t assume only integer values; i.e., the unit 
of time is taken as the stepping interval, and the unit of distance is 
taken as the station separation (center-to-center). Perfect charge 
transfer implies 

q Xil = q x -i, t -i (1) 

and hence unit signal speed. 
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Real devices are typically experimentally characterized by the 
fraction a of the charge which is successfully advanced per step. 4 - 5 If 
a fraction e fails to advance and remains in its original station, then 
the process is described by 

q Xtt =» og.-i.t-i + ■fr.i-i » (2) 

which reduces to equation (1) in the ideal case of a = 1 and « = 0. In 
other words, equation (2) states that the charge in station x at time t, 
q xt , is the successfully transferred fraction of the charge in the previous 
station at the previous time, aq x -i, t -i , plus a fraction e of the charge 
q x ,t-i in station x at the previous time which failed to advance. From 
equation (2) it follows that a fraction £ = 1 — a — « of the charge is 
lost per step. Theoretical expressions for a and/or c are contained ex- 
plicitly or implicitly in the work of several authors on different de- 
vices. 6-9 [In the appendix, equation (2) is generalized slightly to 
include an inhomogeneous term and nonconstant values of a and e.] In 
this context equation (2) was introduced by Berglund to describe in- 
complete transfer in CCD and IGFET bucket-brigade shift registers, 
and he showed an approximate equivalence at low frequencies to an 
(analog) matched transmission line." Traditionally equation (2) is 
associated with probability theory (Bernoulli trials) where a + * = l. 10 
In the case a = c = 1 (outside our primary interpretation) equation 
(2) is usually known as Pascal's triangle although he was preceded by 
Cardan in 1540 who, m turn, cited earlier sources. 11 We give here a 
number of exact and exact limiting (t -> 0) solutions to equation (2) 
useful both in image-detection and shift-register contexts. It should be 
noted that our results can also be applied to Berglund's model 12 of 
bipolar bucket-brigade shift registers which differs from equation (2) 
only through an interchange of the physical interpretation of x and t. 

III. SINUSOIDAL REPRESENTATION 

The dispersion relations can be found from the space and time 
Fourier transforms of equation (2) or by the separation-of-variables 
method. Either approach amounts to seeking a running-wave solution 
of the form 

Ime ,- (u( -*x + „) I 3 | f I i I - 0,1,2, •■• (3) 

where k and w may be complex to account for the attenuation, and f is 
any constant phase factor. Since .r takes on only integer values, equa- 
tion (3) is not affected by adding multiples of 2tt to the real part of k, 
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and thus without loss of generality we require | Re k | ^ tt. In the 
theory of lattice vibrations (where time is continuous but the physical 
system is discrete) this range limitation on Re k is usually expressed by 
saying that the wavelength with one atomic species must be at least 
two lattice spacings. 13 Similarly, since time is discrete, all frequencies 
outside the fundamental range | Re o> | ^ tt are redundant. In the 
theory of sampled-data control systems (where the physical system 
is usually continuous but time is discrete) this is usually expressed 
by saying that the sampling frequency must be at least twice the 
maximum frequency to be detected. 14 (Here the sampling frequency /, 
is once per unit time, i.e., / s = 1 or o>„ = 2tt) . Substituting equation (3) 
into equation (2) yields for the dispersion relation between <o and k 

e iu = ae ik + e. (4) 

Two important special cases of equation (4) are discussed in the rest of 
this section. Since equation (4) shows that the dispersion relations are 
independent of <p, we take <p = hereafter. 

In the case of image detection (or projection), k is real correspond- 
ing to a term in the spatial Fourier representation of the initial image. 
Equation (4) then reduces to 

w (fc) = w ' + «o" = -co*(-fc) (5) 

_! sinfc 

= k — tan j—. — 7- 

cos k + aye 

- i(ln a + \ In [1 + («/a) 2 + 2(e/a) cos k]); (6) 

i.e., the wave is attenuated in time but remains spatially sinusoidal. 

The identity 

_! sin 6 . , -i sin 
tan — ; 1 + tan 



c + cos 6 c~ + cos 9 

was used to express the effect of nonzero e as a separate correction 
term in equation (6). The real and imaginary parts, </ and a>" are 
plotted in Fig. 1. Although not considered here, it is sometimes useful 
to regard a and « as fc-dependent (i.e., wavelength-dependent) quanti- 
ties. In the practical case of small </a, equation (6) is usefully ap- 
proximated by truncating its expansion in powers of e/a after the 
linear term, i.e., 

w(fc) — > k - (e/a) sin k — i[ln a + (e/a) COS fc], e/a — >• 

-> k - e sin k + i[t + e(l - cos k)], a, I -> 0. (7) 
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A phase velocity v' v can be denned which describes the apparent speed 
of the crests and troughs of the attenuating sine wave; i.e., 

vUk, e/a) m ^ - 1 - k- 1 tan" 1 . 8 ' n k . (8) 
* k a/e + cos k 

e sin k . n , . 

-* 1 — r— , «/a-»0. (9) 



a A: 



In particular 



■;(*, e/a) + ■;<*, a/e) = 1 (10) 

»;<*, i) = I (ii) 

y^(fc, e/a) — * 1 as k — > ir if e/a < 1 

-> as k -» 7T if e/a > 1 . 

At long wavelengths (i.e., A: — » 0) co' <-' k while o>" — const ^ /c 2 . Thus 
we ignore the attenuation and find an infinite-wavelength group 
velocity 

^ = »;(o, «/«) - a + «/«r x . 

In the ease of shift-register operation to is real, corresponding to a 
term in the Fourier representation of a time-dependent signal intro- 
duced into one station. Equation (4) then reduces to 

fc( W ) = k' - ik" = -&*(-<■>), g CO ^ 7T (12) 

. . _, sinco 
= co + tan 



— cos CO 

+ t[ln a - $ In (1 + e - 2e cos co)] (13) 

— > co -f- e sin co + i(ln a + e cos co), e — > ,. . 

— > co + e sin co — i[£ + e(l — cos co)], e, I — > 0; 

i.e., the wave is attenuated as a function of x but is purely sinusoidal 
in time at a given x. The identity following equation (6) was again 
used. The quantity ik corresponds to the propagation function of 
(continuous-time) lossy transmission lines; similarly k' and k" corre- 
spond to the phase junction and attenuation junction respectively. 13 
When Fig. 1 is rotated through 180 degrees in its plane, it becomes a 
plot of equation (13). Although not considered here it is sometimes 
useful to regard a and c as co-dependent quantities. 
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Figs, la and lb— For a spatial sine wave of the form exp[i(wi — kx)] with \x\, 
\t\ =6, 1, 2, •• • the complex angular frequency u> = w' + iu>" is plotted versus 
the (real) angular wave number k. After conversion from degrees to radians, k is 
measured in units of (2*- times the stations spacing)- 1 . The curves are param- 
eterized by the ratio e/a where a is the fraction of the charge successfully 
advanced per stepping operation and e is the fraction of the charge which remains 
in a station per step [From equation (6)]. 

Fig. la— The real part of w. (Note that u'(k, e/a) + w'(fc, a/e) = k.) 
Fig. lb— The imaginary part of w. 

After rotating the figures by 180 degrees: For a sinusoidal signal of the form 
expliiut — kxh with \x\, \t\ = 0, 1, 2, • • • (temporal sine wave) the complex 
angular wave number k = k' — ik" is plotted versus the (real) angular signal 
frequency w where w is measured in degrees per stepping interval. The curves 
are parameterized by e where < e < 1 is the physically significant range 
[From equation (13)]. 

Fig. la — The real part of k. 
Fig. lb — The imaginary part of k. 



The phase velocity is given by 



v p (w, «) - fy - 1 1 + w" 1 tan" 1 



sin a? I 
1 — cos wj 



1 - e 



sin co 



(15) 



0. 
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k", and hence the attenuation, is frequency dependent for any c > 
(just as to" is k dependent for any t/<x > 0) ; however, as one would 
expect 16 from causality (Hilbert Transform; Kramers-Kronig Rela- 
tion) in a continuous-time system, the apparent phase velocity w/fc' 
has no dispersionless case for < e < 1 analogous to equation (11). At 
low frequencies k' <— <■> while k" — const — ' u> 2 . Thus we ignore the 
attenuation and find a zero- frequency group velocity 



dk'jO) 
du 



= V,(0, e) = 1 - e. 
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For a numerical example in the small-e approximation, let a time- 
dependent signal be introduced into a shift register at x = and ex- 
tracted at station x. If x is large, the product ex is not necessarily small 
compared to unity; consider the case ex = 0.2. From equation (15) 
the phase transit time through the shift register is x/v,, = xk'/u = 
x + xe sin la/a. Thus the highest frequency component of the signal 
(w = tt) is transmitted in the ideal (i.e., e = 0) time x. Lower frequencies 
take increasingly long with the u = component arriving xe = 0.2 
stepping intervals after the w = if component. The amplitude is at- 
tenuated by the factor e~ k " x = a'e" 0OB a which is a decreasing function 
of frequency. Thus the cut-off frequency (o> = x) is attenuated by 
an additional factor of e~ 2lx = e~ 04 = 0.67 over the attenuation of 
the co = component. 

Equation (4) is invariant under the transformation 

a —> k k —* u) 

a — > 1/a e — > — e/a 

which can therefore be used to derive from each other the two parallel 
sets of relations developed in this section. 

If to is eliminated from equations (3) and (4), equation (25) results; 
if k is eliminated, equation (33) results. 

IV. IMPULSE REPRESENTATION 

Having established the connection with the usual basis for charac- 
terizing continuous linear systems in terms of their effect on sine 
waves, we turn to the impulse representation and investigate the solu- 
tion of a unit charge placed in station x = at t =• with all other 
stations initially uncharged. By considering the spatial distribution 
after the first few steppings as shown in Fig. 2, the solution to equation 
(2) with this boundary condition can be recognized as the binomial ex- 
pansion of (a + c)'; i.e., 

q x . t = G x . t - ('W-, ^ x ^ t 

\x/ (16) 

= otherwise 

where Q is the widely tabulated binomial coefficient. Later we will 
regard G zA as the Green's function for more general initial or boundary 
conditions. Although our interpretation is deterministic rather than 
probabilistic, the terminology and results of probability theory will 
be quite useful. For example, equation (16) is a discrete bivariate 
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Fig. 2 — Evolution of a unit charge initially in station x = at t = 0. A fraction 
a of the charge is advanced per step and a fraction e remains behind. 

(x and I independent discrete variables) distribution, denned by a 
parabolic (partial) finite difference equation [Equation (2)]. 

In the description of practical devices one is usually interested in 
limiting expressions as t and/or e approach zero. Rather than writing 
several limiting forms at each point in the development, we emphasize 
once and for all an approximation which, in one variation or another, 
is frequently useful, 



' = n _ /V = *< m <>-<> 



(a + e)' = (1 - I) 1 = e 



tf « 1. 



(17) 



Although in a real device the loss per transfer ( may be orders of mag- 
nitude less than unity, the number of stations and hence transfers 
may be so large that tt is comparable to unity, and the further ap- 
proximation e~" f& 1 — tC is not generally accurate. 

Consider first the spatial distribution of charge with time regarded 
merely as a parameter, i.e., a horizontal row in Fig. 2. In probability 
theory it is usually assumed that a + e = 1 (probability of heads 
plus probability of tails are unity). To make direct use of extensively 
compiled properties define 



a' = «/(l - 0, *' - «/(l - 



(18) 
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so that a' + e' = 1. Then equation (16) becomes 

G Xll = (1 - ^'(^)a'V'- 1 = (1 - tyb z , x - 0, 1, ■ • • , t (19) 
which defines b x the binomial distribution 17 with normalization 

Mo - £, b z = (a' + e')' - 1 (20) 

i-0 

and (1 - l) 1 is just a position-independent scale factor. (The inter- 
mediate steps involved in carrying out the numerous sums which 
follow can be found in many texts on probability theory. 1 ) 

If the speed of the pulse is defined as the speed of the center of 
mass of that charge remaining at any time, then from 

(x) m j^ **>* = «'< < 21) 

i = 

it follows that the speed is a' = (1 + c/a)" 1 , i.e., the same as the 
infinite-wavelength group velocity following equation (11). In general, 
the pulse can be characterized by its central moments m of which the 
rth 

Mr = E(s- {x)yb It r = 0, 1, ••• , co (22) 

x = 

can be obtained explicitly as the coefficient of ff/r\ in the Taylor- 
series expansion in s of the central-moment- generating junction 11 

e- {z) J2e' x b x = e- (I> (e'+aV)\ (23) 

i = 

Thus in particular the variance of the pulse, a useful measure of its 
spread, is fi 2 = tut t' = t(2 + a/e -I- e/a) -1 where the last form again 
explicitly exhibits the fact that normalization-independent properties 
depend on a and i only through their ratio. Higher central moments 
are conveniently obtained from Romanovsky's recursion relation. 17 - 18 

Mr+1 = aV(frM,-, + dnr/da% (24) 

Some further useful relations are as follows: When the pulse spreads 
excessively for digital applications, one alternative is to utilize 
2p + 1 consecutive stations for one pulse. It then follows directly 
from the parallel-axis theorem for the moment of inertia that the 
variance is increased by (only) the additive term p(p + l)/3 over 
the previous result. 
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At time t the ratio of the charge in the next-to-leading station to 
that in the leading station of the pulse is UJa. Thus the leading station 
of the pulse is also the most heavily charged until time t = a/e when 
the peak starts to drop back from the lead. At all times the charge 
distribution falls off monotonically from the peak. At any time the 
fraction of the remaining charge confined to the leading station is 

(1 + e/a)-K 

The results in Section III can be derived as consequences of those 
of Section IV, and conversely. For example the single-spatial-Fourier- 
transform solution of equation (2) can be obtained via the binomial 
distribution regarded as a Green's function (in the terminology of 
mathematical physics) . Thus the initial sinusoidal spatial charge dis- 
tribution q Xj0 = Im er tk * evolves into q Xit at time t where 

g,., = Im T,G z , A e- ik{z - z,) 

(25) 
= Im (e + ae ik )'e- ikl = Im (a + «f ") V*" - *' . 

Equation (25) is equivalent to equations (3) and (4) with <•> elimi- 
nated. Equation (6) results if the real and imaginary parts of equation 
(25) are extracted. Apart from the normalization, the evolution factor 
which carries the initial spatial distribution e~^ hx into the later distribu- 
tion can be recognized as W + a e ,-fc )', the characteristic junction 
(Fourier transform) of the binomial distribution in probability 
theory. 17 In the last form of equation (25), (a + ^e~ ik ) t evolves the 
unperturbed or ideal (« = 0, a = 1) solution e ik{t ~ x) into the actual 
solution. Finally (1 + (c/a) e _fA ) ' (— » exp {e^He/a) as«/a-»0) evolves 
the e = solution a*e**<^*>. 

As before [cf. equation (7)] we seek an approximation of the Green's 
function in the limit that e/a, but not necessarily te/a, is small compared 
to unity. If I is also small, repeated use of equation (17) yields 

G Xit — * e~"p(t — x), e, I — * 0, t — > oo t tt, tl = const. 

\ = te, x = t, t - 1, ••■ (26) 

where p{u) = e -x A u /w! (w = 0, 1, • • •) is the Poisson distribution. 17 
Sum rules similar to those illustrated for the binomial distribution 
can be carried out yielding in particular a pulse speed (x)/t = 1 — 
\Jt = 1 — e and a variance A = U. As long as the pulse peaks near 
or at the leading station, the Poisson distribution is usefully accurate 
for the more strongly charged stations even when t is only modestly 
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large. A comparison of the binomial distribution and its Poisson ap- 
proximation are given in Table I. 

It is impractical to summarize the very many relevant exact and 
approximate properties of the binomial distribution here, particularly 
since a thorough compilation is already available." We have, however, 
given a few of the most important properties to illustrate the useful- 
ness of this connection with the extensive literature of probability 
theory. 

Turning now to a shift-register interpretation of equation (16), we 
regard z as a parameter (the number of the stations which are used as 
the register) and study the time-dependence of the charge in station 
X — 1, i.e., a vertical column in Fig. 2. The development proceeds in 
close analogy to that of the binomial distribution and most of the 
motivating remarks need not be repeated. 



Table I — A Comparison of the Binomial and Negative Binomial 
Distributions with their Poisson Approximation 
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For the case a = 0.9 (fraction of charge successfully advanced per step), e =0.1 
(fraction which remains in its station per step), and I = (fraction of charge lost 
per step). 

Column 1 : Station number x for columns 2, 3, and 4. 

Column 2: Initial distribution of charge among the stations. 

Column 3 : Binomial distribution of charge among stations eight steps later If rom 
equation (16) and (19)]. . 

Column 6: Detection time t for Columns 5 and 4 as measured from time ot initial 
detection, t = 0. 

Column 5: a times the negative binomial distribution of charge nondestructively 
observed in station 7 (a times equation (27)) as a result of initial spatial 
distribution of Column 2 = distribution of charge observed by a charge- 
removing detector in station 8 (cf. Appendix). 

Column 4: Poisson approximation to Columns 3 and 5 [from equation (26) or (34)J. 
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It is notationally convenient to use a new time t which starts at 
time t = x — I when some of the charge placed initially in station 
zero is first observed in station x — 1, i.e., t = t — x + 1. Then equation 
(16) becomes 

GU., = ( T ^*f 7 ^"V, r = 0, 1, ••• ,00. (27) 

If a T is defined as 1 — « (so that a + + e = 1), equation (27) becomes 
_,/ a V/r + X - l\ t, T 



*—-«■' Ir^iA ;ii> v — ^t±-> < 28 > 

and defines ?? T the negative binomial distribution 17 with normalization 
Mo = f> r = a tx (l - e)" = 1. (29) 

T = 

Since charge which fails to advance is observed again, the sum of all 
charge (nondestructively) observed in station x — 1, i.e., a I_1 (l — e)~", 
may exceed unity even if I ^ (cf. the discussion of boundary con- 
ditions in the appendix). Whereas all normalization-independent 
properties of the binomial distribution depend on a and e only through 
their ratio e/a, here all normalization-independent properties depend 
only on e (cf. Figs. 1 and 2). 

The charge in station x — 1 is observed at a mean time (t) = x — 
1 + (t) where 

< T ) = f; m r = . C€ (i _ e )-' (30) 

and thus in the shift-register context the pulse speed can be defined 
as (x — \)/{t) which approaches 1 — e for large x; i.e., the same as 
the zero-frequency group velocity following equation (15) but faster 
than the definition following equations (11) and (21) except in the 
case of no charge loss (f = 0) when all definitions agree. 

The rth central moment of the charge sequence in station x — 1 is 

M, - Z (r - (r)Yn T (31) 

T = 

which can be obtained by direct calculation or as the coefficient of 
s r /r\ in the Taylor expansion in powers of s of the central-moment 
generating function 17 
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e" iT) E e T n r = <T* <T, (1 - (*T(1 - *Y- ( 32 > 

T = 

In particular, the variance is ^ = xe (1 — <)~ 2 . 

The first charge observed in station x — 1 is a fraction (1 - c) a of 
all that will be observed in that station. It is also the largest charge 
observed in x — 1 if Xe, the ratio of the second amount observed to 
the first, is less than unity. Thus the input, Xe, for the example at the 
end of Section III can be obtained experimentally from the ratio of 
the first two nonzero charges present in x - 1 due to a single initial 
charge in x = 0. More accurate graphical methods actually developed 
for statistical contexts can be directly applied here to infer the param- 
eters from the final distribution resulting from an initial pulse." 

According to equation (2) a sinusoidal signal q , t = Im e iwt present 
in station zero at time t adds an additional charge &qi,t+i to station 1 
at t + 1 where Aqi,t+i = «9o,* (cf. the Appendix) . Thus from equation 
(28) the charge present in station x at time t is 



t T~* n iuU-z-T) 

<lx, t = Im 2-/ (j x -i.T<xe 



= Iml-^— )«'-' = Im' ^ l 



(33) 



4 - «f 

which is equivalent to equations (3) and (4) with k eliminated. Equa- 
tion (13) results if the real and imaginary parts of equation (33) are 
extracted. Apart from the normalization the transfer factor which 
multiplies the initial signal e iwt in equation (33) is the characteristic 
function of the negative binomial distribution. 

In the limit of small t and t, but not necessarily small ex or d, G XiT 
becomes exp (-tz) times the Poisson distribution; i.e., 

g mr _» e - V a A t /t!, t, e -* 0; x -* « ; te, *x = const A = xe (34) 

which is functionally equivalent to equation (26) . Thus as illustrated 
in Table I, the distinction between the spatial and temporal projec- 
tions of the impulse solution disappears in the Poisson limit [cf. equa- 
tions (7) and (14)]. 

V. SUMMARY 

Without reviewing any one of the three specialized fields individ- 
ually, we have emphasized by way of an explicit case of current in- 
terest the connection between, on the one hand, the present practice 
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of characterizing discrete charge-transfer devices in terms of their 
charge-transfer efficiencies and, on the other hand, the two well- 
developed fields of (analog) sinusoidal signal analysis and probability 
theory (interpreted deterministically) . 

Parenthetically we note that the model of this paper [equation (2)] 
is of some tutorial interest in that it permits the basic ideas of a 
traveling-wave description of discrete-space-and-time linear systems 
(including therefore as special cases the limits of continuous space 
and time) to be exemplified immediately in a unified manner via 
a simple device. The usual textbook vehicles such as sampled-data 
control systems, lattices, or transmission lines suffer tutorially in 
varying degrees from being too restricted (discrete in only space or 
time; no attenuation) and requiring a lengthy and specialized physical 
explanation of the origin of the basic equation in a less easily 
visualized system. 
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APPENDIX 

In the main portion of this paper the boundary conditions were 
deemphasized. For example, in shift-register operation the charge 
distribution was assumed to be unaffected by its own detection. Let 
q x , t be the nondestructively observed charge as distinguished from 
q' zt , the charge measured by a destruction process that removes all 
of the charge from station x. Then by equation (2) these two limiting 
cases are simply related as follows 

Qx.t = ag.-i.i-i ; (35) 

i.e., the expressions given for a nondestructively observed shift- 
register need only be multiplied by a to obtain their destructively 
measured counterparts in a register with one more station. In the 
case of a detection process which removes a fixed fraction of the 
charge, the signal can be obtained from the difference between suc- 
cessive measurements. 

Similarly the two limiting cases at the input station of the register 
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are Dirichlet (voltage) and Neumann (charge) boundary conditions. 
In the Neumann case a fixed amount of charge is added to whatever 
charge may already be in the station. In the Dirichlet case a fixed 
applied voltage adds charge as necessary to what is already in the 
station thereby achieving a total charge predetermined by the voltage 
level. By equation (2), and in analogy to equation (35), these two 
cases are simply related because a fixed total charge q x , t in x at t 
injects a fixed additional charge aq Xit into x + 1 at t + 1 (cf. equa- 
tion (33)). 

To generalize equation (2) consider a constant inhomogeneous term 
q d which might, for example, represent the amount of dark-current 
charge added per station per stepping interval. Then (2) becomes 

q Zit = aqz-i.t-! + 6&,,-i + q d ■ (36) 

An important particular solution of equation (36) is the spatially- 
uniform rising-in-time solution which starts from a spatially uniform 
distribution q at t = 0: 

1 — 0' 
q z .t = ff'go + x _ fl g" /3 = a+e, t = 0, 1, • • ■ (3?) 

= q + q*t, 1=0. 
This could describe a recirculating memory where the output of the 
last station is fed back into the initial station. The constant-in-time 
rising-in-space solution is 

9"-**+iL.\-l *, Y-rh' '-I 1 *- (38) 

= 9° + T^l qd ' * = °' 

where q is now the charge in the (initial) station x = 0. This steady- 
state distribution could be maintained in a shift register by removing 
a charge q d — (1 — e)q from the initial station per step. Equation (2) 
is now understood to describe a signal (homogeneous solution) imposed 
on' some particular solution to equation (36) describing the back- 
ground charge. 

Next assume, as is certainly true in some devices, that the transfer 
fractions a and e depend upon at least some of the g's; i.e., the process 
is nonlinear. Then in the small-signal case (whether or not q d is fur- 
ther generalized to be a function of x and t) a and e in equation (2) 
would still be independent of the magnitude, of the signal charge but 
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no longer independent of x and t even though all stations were phys- 
ically equivalent apart from their momentary charges. In the case of 
equation (37) a and e would still be x independent; in the case of 
equation (38) still t independent 10 (cf. nonuniform transmission lines). 

Before proceeding we switch to the suitably developed terminology 
of probability theory (without necessarily implying a probabilistic 
interpretation). Actually, this has some plausible justification. To the 
extent that equation (2) applies independently to each of the indi- 
visible electrons or holes in the charge, a and e must be interpreted as 
transition probabilities, and our deterministic interpretation results 
only because the large number of electrons or holes permits fluctuations 
to be ignored. In this probabilistic sense one can define the information 
content of a discrete charge distribution. 20 Similarly in equation (36) , 
q a (x, t) could be a random variable describing the introduction of 
extraneous noise. 

To complete the small-signal or linear model for the signal we note 
that in practice the stations are driven by w-phase time-dependent 
voltages such that in a coordinate system y which moves at the ideal 
signal speed (y = x — t) the voltage distribution among the stations 
appears constant in discrete time and periodic in discrete space with 
a period of n stations. In each cycle of n stations one station is the 
designated potential well which, in the ideal case, carries all of the 
charge in that spatial cycle. The distortion of the signal corresponds 
to the transitions which the electrons or holes make to neighboring 
stations and, eventually, to neighboring wells. In a causal system the 
homogeneous equation for the signal charge takes the form 



Qy,> 



= Z E T u , v . ,,,rq v ',r (39) 



where the transition elements of the T matrix are g-independent in a 
small-signal theory. (In some cases, the upper limit of 3 is t rather 
than t — 1.) 

If the system is memoryless such that the spatial charge distribution 
at one time is sufficient to determine the distribution at the next time 
and hence for all time, then equation (39) reduces to the Markoffian 
form 

Qy.l = 2 ?».»'. f.t-lfr'.t-l • (40) 

v' = — oo 

This simplification corresponds, for example, to neglecting traps which 
accept signal charge at a rate proportional to the local signal and then 
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emit the charges into the signal over a time scale long compared to 
the stepping interval. 

When the background charge distribution, if any, is time independent, 
T u ,y t ,r depends on t and t' only through their difference. If the back- 
ground distribution is spatially constant, T is a periodic function of 
y _ y'- } i e-) t describes a random-walk" 1 problem on a one-dimensional 
lattice with anisotropic and spatially periodic transition probabilities 
(cf. Floquet's theorem 13 ). 

Finally, of course, large signals under nonlinear conditions are of 
interest. However, it is not yet clear that there is any dependence of T 
on the q's which is more general than the specific device where it 
arises. 

In the moving coordinate system equation (2) becomes 

Qv.t = <xq u .t-i + «9i/+i.<-i ( 41 ) 

which is functionally equivalent to equation (2) and distinct only in 
that the roles of a and «, as well as the sense of spatial direction, are 
reversed. If the transition probabilities (T matrix elements) are small 
out of the wells and large out of the other stations between wells, then 
this 1-phase equation can be used to approximate an n-phase system. 
The unit of time is taken as the effective stepping interval equal to n 
actual stepping intervals, and successive integer values of x or y label 
the n-cycle groups or the wells. The matrix elements e and a in equa- 
tion (41) then describe a net transfer of charge between, or retention of 
charge by, the n-cycle moving groups during one (effective) time unit. 
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